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(Borexino Gollaboration) 

Borexino is a liquid scintillation detector located deep underground at the Laboratori Nazionali 
del Gran Sasso (LNGS, Italy). Thanks to the unmatched radio-purity of the scintillator, and to the 
well understood detector response at low energy, a new limit on the stability of the electron for decay 
into a neutrino and a single mono-energetic photon was obtained. This new bound, r > 6.6 x yr 
at 90 % G.L., is two orders of magnitude better than the previous limit. 


The conservation of electric charge, suggested since 
the 19th century, is fundamental to the physics of the 
Standard Model as a direct consequence of Maxwell’s 
equations and the unbroken U{1) gauge symmetry of the 
electro-weak theory. Despite the present undisputed va¬ 


lidity of this law, experimental tests of charge conserva¬ 
tion remain a way to search for physics beyond the Stan¬ 
dard Model, and deserve to be investigated with the high¬ 
est possible sensitivity. An experimental search for the 
hypothetical charge non-conserving decay of the electron. 
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which is the lightest known charged particle, into a neu¬ 
trino and a photon is reported in this paper. No presently 
viable theory predicts such a decay, and a large charge vi¬ 
olation is excluded by the absence of macroscopic effects 
in matter. The electron decay is, however, discussed in 
the literature, e.g. [H-Q and references therein. 

The Borexino detector Q is a unique tool to undertake 
a search for electron decay. The unmatched radio-purity 
of the liquid scintillator (a mixture of PC and PPO Q) 
was obtained by means of a very careful selection and 
cleaning of materials, the development of innovative pu¬ 
rification techniques , and the extreme care taken 

during procurement and handling of the scintillator fluid 
and the filling of the detector. In addition, the successful 
measurement of several solar neutrino components with 
Borexino (^Be 11. , [ 13 , ®B [ 13 , pep [ 3 , and pp M) 
demonstrates an unprecedented sensitivity which can be 
applied to the search for a two-body e —>■ + 7 decay 
resulting in the emission of a 256 keV photon. 

The Borexino collaboration has already searched for 
this decay using the Counting Test Facility (CTF) 0 , 
a 4 m^ liquid scintillator prototype detector operated to 
prove the feasibility of Borexino and to develop the nec¬ 
essary purification techniques for the scintillator. This 
earlier measurement used PXE scintillator 17| , and set a 
lower limit for the lifetime of this electron decay mode of 
4.6 X 10^® yr [l3. The work reported here uses a mixture 
of PC and PPO scintillating solution in the Borexino de¬ 
tector instead. A sensitivity improvement of more than 
two orders of magnitude is obtained thanks to a larger 
fiducial mass, higher statistics, lower background, and 
significantly improved data analysis. The latter is the re¬ 
sult of a deeper understanding of the detector energy re¬ 
sponse, obtained by extensive calibration campaigns [l 3 
and encoded in a tuned Geant4 simulation. 

In Borexino, neutrino interactions and background 
events are detected by means of scintillation light col¬ 
lected by some 2212 photomultipliers (PMTs) (various 
PMTs were taken offline over the years of Borexino opera¬ 
tions due to electronic failure). The PMTs are supported 
by a Stainless Steel Sphere (SSS). The sphere contains, a 
thin concentric nylon vessel which separates the detection 
scintillator fluid, a PC-|-PPO(1.5 g/1) scintillating mix¬ 
ture, from an outer buffer liquid, a mixture of PC and 
DMP without PPO. A second larger nylon vessel acts as 
a Radon barrier. The SSS is contained in a large Water 
Tank (WT) instrumented with PMTs, shielding the in¬ 
ner spheres against external radiation and providing an 
active muon Cherenkov veto. For more details about the 
Borexino detector see (3 and 

The light collected by all the PMTs provides four ba¬ 
sic measurements: 1 ) the total deposited energy, recon¬ 
structed from the number of collected photo-electrons; 
2 ) the event position, reconstructed by means of a time- 
of-flight fit between the PMT signals; 3) the particle 
type (a-like or /3-like or muon), reconstructed from pulse 


shape analysis; and 4) the WT Cherenkov signal. For 
more details see 121 . 

The analysis performed for this paper has taken advan¬ 
tage of the recent measurement of pp solar neutrinos 0 , 
as the expected photon visible energy from the hypothet¬ 
ical two-body electron decay occurs in the region where 
the pp neutrino signal is dominant. The correlation be¬ 
tween pp neutrinos and the electron decay signal is dis¬ 
cussed below. The study of the lower energy region of the 
spectrum (165-590 keV) poses special challenges related 
to the linearity of the energy response, light quenching 
processes, and ^“^C event pile-up. These were tackled for 
the measurement of pp neutrinos for which special tools 
were developed, which are also applied in this work. 

Data used in this search were collected during the 
“Phase 2” operations of the Borexino experiment, started 
in 2012 after a set of calibration runs and more than a 
year of liquid scintillator purification. The calibration 
was performed with internal 7 , /3, a, and neutron sources, 
which yielded a meticulous understanding of the detector 
response over a large energy range. Scintillator purifi¬ 
cation, done by means of water extraction and nitrogen 
stripping, substantially reduced radioactive backgrounds. 
In particular, ®®Kr concentration is now compatible with 
zero (from ^35 cpd/100 t in Phase 1), and the ^^®Bi con¬ 
tent was reduced by a factor ^4, to about 20 cpd /100 
t. 232rp]^ concentrations were at a record low, 

< 10“^® g/g. Data for this work were acquired from Jan¬ 
uary 2012 to May 2013, corresponding to 408 live days. 

The energy spectrum used for the electron decay search 
is shown in Fig. [T] (black points with error bars) to¬ 
gether with the main fitted components (color lines). At 
low energy (below 200 keV) the count rate is dominated 
by the /3-decays of ^^C, with a measured abundance of 
(2.7 ± 0.1) X 10“^® g/g [I 3 with respect to ^^C. The 
mono-energetic peak in the central part of the spectrum 
corresponds to 5.3 MeV a—particles from ^^®Po decay, 
which shifts downward to approximately 400 keV elec¬ 
tron equivalent by quenching in the scintillator 2]| . The 
arrow indicates the position of the hypothetical 256 keV 
7 peak from electron decay. 

The event rate from ^"‘C, while intrinsically low, still 
yields large event statistics over the entire scintillator 
volume, on the order of 5 x 10® events for the lowest 
energy bin used in this analysis. This requires very pre¬ 
cise fitting models to keep the systematic uncertainties 
at or below the statistical fluctuations of approximately 
0.14%. A consequence of ^"^C decays is a non-negligible 
occurrence of pile-up events, when two or more indepen¬ 
dent decays (mostly ^^C) occur close enough in time not 
to be separated. Two events can be distinguished with 
>50% efficiency when they are separated by more than 
230 ns. The energy spectrum of pile-up events in the re¬ 
gion of interest above the ^'’’C end-point is similar to the 
electron recoil spectrum induced by pp solar neutrinos. 
The pile-up spectrum obtained for the measurement of 
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pp neutrinos [l5l| . shown in magenta in Figll] is included 
as a separate component in the spectral fit used in this 
analysis and labeled “synthetic” pile-up. The synthetic 
pile-up spectrum is constructed by overlapping triggered 
events with PMT hits recorded in the tail of the acqui¬ 
sition gate of such events, well after the triggered scin¬ 
tillation pulse has decayed away. These late PMT hits 
represent a data-driven, threshold-less and random sam¬ 
ple of activity (data -I- noise) in the detector. Added to 
triggered events they boost the pile-up contribution by 
a known amount, allowing the pile-up spectrum to be 
extracted. 

The scintillation signal from a 256 keV 7 produced in 
the e~ ^ + Ve decay is equivalent to one produced by a 
220 ±0.4 keV electron [l^. The energy shift is due to the 
partial light loss from quenching, i.e. the non-linearity of 
the scintillation response with electron energy. Quench¬ 
ing is modeled by the standard Birks’ formalism 2l|, 
which relates the density of light production ^ to the 
ionisation density 


dL 

dx 


dE 

dx 




( 1 ) 


where k and B are the Birks’ parameters. The average 
number of collected photo-electrons (p.e.) Q produced 
by an electron of energy E can be obtained by integrating 
O- It is convenient to present the result in the form: 


Q = LY-E-fikB,E), (2) 


where LY is the light yield for electrons expressed in 
p.e.-MeV“^ and f{kB,E) is a light deficit function, i.e. 
the result of the integration of eq. [T] along the path nor¬ 
malised to unity at 1 MeV. The f{kB,E) is a mono- 
tonically increasing function in the region of interest. 
Thus, the average light yield of a 7 absorbed by multiple 
Compton scatterings at low energy followed by photo¬ 
absorption, is lower than that released by a single elec¬ 
tron of the same initial energy. This fact is crucial be¬ 
cause the quenched 256 keV -7 energy partially overlaps 
with the tail, which then required a special analysis. 

The number n of PMTs that give a valid hit within a 
time window of 230 ns is approximately proportional to 
the energy deposit and therefore to the total charge Q col¬ 
lected by the PMTs. The relation between n and Q is: 


Npmt log(l- 
l+5-log(l-lv^) ’ 


( 3 ) 


where Npmt is the total number of PMTs of the detec¬ 
tor, and Pc is a geometric correction factor obtained by 
means of MC. This approach is the same used in solar 
neutrino analysis and is described in detail in 1 ^ and 
The statistics of n is easier to model with respect 
to the number of collected photoelectrons Q because the 
statistical distribution of multiple hits in PMTs depends 


Fit result for the electron decay rate = 1.23 cpd/100 tons 



FIG. 1. Energy spectrum between 150 and 600 keV. The most 
prominent features are the /3-spectrum (green), the peak 
at about 400 keV from ^^°Po o-decays, and the solar neutri¬ 
nos, grouped in the blue curve except the crucial pp neutrinos, 
which are shown in cyan. The effect of event pile-up, mostly 
overlapping events, is shown in dashed-pink. The hypo¬ 
thetical mono-energetic 256 keV 7 line is shown in red at its 
90% exclusion C.L. with an arrow indicating the mean value 
of the detected energy, which is lower than 256 keV because 
of quenching. The fit is done in the range 164-590 keV. 


on the details of the electronic response, and it is not 
known with sufficient precision. 

The most crucial part of the analysis is the behaviour 
of the energy resolution as a function of the energy. Thus 
the variance of the energy resolution (in terms of the used 
energy estimator) is modeled as [ 2 ^ : 

= N{po-piVi) + n^{vT{n) + Vf{N)) + al + a^^t (4) 

where po = 1 — pi, N = {f{t))rp is the average number 
of operating PMTs during the period of the data acqui¬ 
sition, and T, f{t) is the number of operating PMTs 
as a function of time, normalised as /(O) = Nq. Here 
Nq is the number of working PMTs at the beginning, 
Vf{N) = (/^(t))^ — {f{t))‘^ is the variance of f{t) over 
the period of data acquisition, <Jd is the contribution of 
the dark noise (fixed at the measured value as an average 
over all PMTs) and, dint is the contribution due to the 
smearing of the intrinsic line. 

The probability pi that a single PMT is hit in a given 
event depends on event position, and is related to the 
energy estimator, n, as n = N ■ pi. The parameter, vi = 
^ (Pi(^)'^i('^))yi is the weighted average variance of 
the single PMT response, where ui(l^) is defined as the 
variance of the pi (averaged over all PMTs) for an event 
whose position is 7^. The quantity, ui, is calculated over 
the detector’s Fiducial Volume (FV) using Monte Carlo 
(MC) simulations. Due to the very narrow energy region 
of interest, - essentially the region of the ^^C tail -, the 
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energy dependence of vi in this analysis can be neglected. 
From the above the value of = 0.17 is computed. 

Finally, vt is the variance of the number of triggered 
PMTs over the FV for a fixed energy. It takes into ac¬ 
count the non-uniform light collection over the detector’s 
volume and the additional variability of the number of 
triggered PMTs in different locations within the FV. It 
is energy-dependent and was found to be proportional to 
energy in the region of interest. Its energy dependence 
was modelled as vt = v^n [ 13 , leaving the constant of 
proportionality as a free parameter of the fit. 

The shape of the scintillation line (i.e., the energy re¬ 
sponse of the detector to a mono-energetic source uni¬ 
formly distributed within FV) is another important com¬ 
ponent of this analysis. The familiar gaussian approxi¬ 
mation fails to describe the tails of MC-generated mono- 
energetic peaks, even for the statistics on the order of 
10^ events. In the previous ^Be solar neutrino analysis 


llj this problem was solved by using a generalised F- 
function 2^ to fit the mono-energetic ^^'^Po peak. How¬ 
ever, while the quality of the fit to the ^^°Po peak is 
insensitive to the residual deviations in the tails, this is 
not the case for the spectrum as all the events in the 
fraction of the spectrum above the end-point 
originate from spectral smearing so that the statistics in 
the tails are much higher. For these reasons, a different 
procedure was adopted. 

The ideal detector response to a point-like mono- 
energetic source in the center is an exact binomial dis¬ 
tribution which can be well approximated by a Poisson 
distribution. However, the ideal width of the distribution 
must be modihed to include the additional spreading of 
the signal due to various factors. The problem, both 
with binomial distribution and with its Poisson approxi¬ 
mation, is that their width is defined by the mean value, 
which introduces an unwanted and unphysical correlation 
between the position of the peak and its width. To attack 
this problem, the response function was approximated by 
a Scaled Poisson Distribution (SPD) defined by: 


sx 

fix) = (5) 

(sx)! 

where x is the independent variable whose mean value is 
n and variance cr„. This function has two free parameters 
/i and s, which can be evaluated using an expected mean 
and variance. The agreement of this approximation with 
the detector response function was tested with the Borex- 
ino MC model. It was found at low energies that the 
function (l5|) reproduces the scintillation line shape much 
better than a generalised F-function (GGF) up to statis¬ 
tics of 10 ® events per bin, while at energies just above 
the tail both distributions give comparable results. 
Therefore the SPD approximation was adopted, and the 
quality of the fit was estimated using a criterion. As 
an example, with 10^ mono-energetic events for n=50 
(approximately 140 keV) we found x^/n.d./.=88.0/61 for 


the GGF compared to ^/n.d.f.=b9.‘i/Ql for the SPD. 
In this example the events were uniformly distributed in 
the detector before the FV is selected. 

As proven by MC calculations, the SPD works well in 
the region of interest despite the additional smearing due 
to the aforementioned factors. This is a result of the fold¬ 
ing of the relatively narrow non-statistical distributions 
by the much wider base function. The MC shows that 
such an absorption results in the smearing of the total 
distribution without changing its shape. 

Only a fraction of the total response for the mono- 
energetic 256 keV 7 enters into the analysis window above 
threshold, which makes the signal look similar to the 
pp-spectrum and produces a strong correlation between 
them (see Fig. 1). Note the correlation between electron 
decay and pp solar neutrinos: with an unconstrained pp 
flux, 12 cpd /100 t of electron decay candidates corre¬ 
spond to about 134 cpd/100 t from pp solar neutrinos in 
the fit. To break the degeneracy, the pp-neutrino rate can 
be constrained either by the value measured by experi¬ 
ments other than Borexino or at that predicted by LMA- 
MSW theory. We chose to use data from radiochemical 
experiments only to obtain a mo del-independent result. 

The deviation from the observed pp-neutrinos rate is 
controlled by a penalty term added to the in the form: 

2 ^ {pp - {PP)f 

^PP rr2 


where pp is the pp-neutrino count rate found by the fit, 
(pp)=134 cpd/100 t is the rate expected in Borexino, 
and (Tpp=13.3 cpd/100 t its variance. These are calcu¬ 
lated using the solar pp-neutrino flux obtained from the 
combined analysis of the GALLEX/GNO and SAGE ex¬ 
periments 2^. The parameter (Tp’y=2.7 takes into ac¬ 
count the systematic error (2% as reported in [l5|) due 
to the uncertainty in the FV mass. 

A search for the 256 keV 7 line was then undertaken 
using “standard” fit conditions, dehned by: a) an en¬ 
ergy estimator: number of triggered PMTs in a fixed 
time window of 230 ns (npmts); b) a fit range: 62- 
220 npmts, corresponding to 164-590 keV; c) a fiducial 
mass: 75.5 t (R<3.02 m and \Z\ <1.67 m) which is the 
same as in solar ^Be analysis. The values of solar neu¬ 
trino rates other than pp are constrained either at the 
results found by Borexino in a different energy range 
of R(^Be)=48±2.3 cpd/100 t [llj, or fixed at the pre¬ 
dictions of the SSM in the MSW/LMA oscillation sce¬ 
nario, R(pep)= 2.80 cpd/100 1 , R(GNO)=5.36 cpd/100 1 . 
The rate was constrained at the value found 

in an independent measurement R(^^C)=40±1 Bq (or 
R(i4c)=(3.456±0.0864)xl0® cpd/100 t). The pile-up 
rate was constrained at the values found with the above 
algorithm to be 321 ± 7 cpd/100 t. Other background 
components were left free (®®Kr,^^°Bi and ^^°Po) and the 
rate of ^^^Pb R(^^^Pb)=0.06 cpd/100 t was calculated by 
means of identified ^^^Rn events, which were measured 
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by the detection of coincidences. The light 

yield and the two energy resolution parameters {vt and 
aint) are left free in the fit. The position of the ^^°Po 
peak is also left free, and decoupled from the energy scale 
because of its poorly known quenching factor. 

The detection efficiency of 256 keV 7 s was determined 
by means of the MC code. A set of uniformly distributed 
7 s were simulated inside the entire IV. The events pass¬ 
ing the same set of cuts used for real data selection and 
reconstructed within the FV were used to determine a 
global efficiency of e = 0.264, which includes FV cut. 

The upper limit S for the given confidence level C.L. 
comes from the following integral: 

pS poo 

/ /(e)d^ = C.L. / /(Od^, ( 6 ) 

Jo Jo 

where /(^) is the probability distribution function. With 
this procedure, we found the 90% C.L. upper limit S=379 
events, in a period T=408 d and fiducial mass of 75.5 t. 

The corresponding lifetime tbx > 7.2 x 10^® yr (sta¬ 
tistical only) was obtained using the relation: 


where Ne = 9.19 x 10®^ are the electrons in the IV mass 
of 278 t and e=0.264 is the detection efficiency. 

The systematic error depends on three factors: the 
choice of the energy estimator, the value of the quenching 
parameters, and the knowledge of the fiducial volume. 

The main one is the choice of the energy estimator, and 
in particular the length of the 230 ns time window, which 
affects the relevance of pilc-up events and dark noise hits. 
We have throughly studied this effect by means of MC 
simulations. It amounts, for this analysis, to about 8 %. 

The second one is the uncertainty is the position of 
the 256 keV 7 -peak with respect to the ^^C spectrum 
depends on the quenching factor, kB, and the pre¬ 
cise position of the end-point of the ^^C /3-spectrum. 
Since the underlying ^"‘C /3-spectrum is falling expo¬ 
nentially, a shift of the peak position to lower energy 
decreases the sensitivity of the search and vice versa. 
In the present analysis the kB parameter is fixed at 
kB = 0.0109 cm • MeV“^, determined from calibration 
data with mono-energetic 7 sources. In order to study 
the related systematics in a model independent way, the 
position of the 7 -peak within the bounds was obtained 
by calibration (1%) instead of varying the kB value. It 
should be noted that no degradation of the fit quality 
for a wide range of kB values was observed, due to the 
absence of known 7 -sources in the region of interest and, 
therefore, in the fit model. The final result is that kB 
error adds less than 2 % to the total systematic error. 

The third systematic error is the knowledge of the FV 
mass which, at the energies of interest, adds an uncer¬ 
tainty of less than 1 % on electron decay rate. 


The three effects described above are correlated. The 
total systematic error was obtained by building a set of 
probability profiles for a corresponding set of input pa¬ 
rameters. The probability profiles were renormalised to 
physical regions, excluding non-physical values in case 
of a negative number of candidates where the region be¬ 
low zero was cut and the p.d.f. was renormalised to unity. 
The final probability profile was obtained as the weighted 
sum of the separate p.d.f.s, with the weights being the 
probability of occurrence of the corresponding value of 
the parameter. The normal distribution was used for 
both the error on the FV and for the relative shift of the 
256 keV 7 with respect to the electrons scale. The ad¬ 
dition of the total systematic error degrades the purely 
statistical result by about 8 %. 

The lower limit on the mean lifetime for the decay 
e~ —7 -I- z/ is Tbx > 6.6 x 10^® yr (90% C.L.), which 
improves the existing limit from CTF 0 by more than 
two orders of magnitude. The sensitivity is such that a 
5(7 discovery signal would have been possibile with an 
electron lifetime of 1.9 x 10^® yr. 
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